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Abstract 

The inclusive Higgs production cross section from gluon fusion is calculated through 
NNLO QCD, including its top quark mass dependence. This is achieved through a 
matching of the 1/M t expansion of the partonic cross sections to the exact large- s 
limits which are derived from fc^-factorization. The accuracy of this procedure is esti- 
mated to be better than 1% for the hadronic cross section. The final result is shown to 
be within 1% of the commonly used effective theory approach, thus confirming earlier 
findings. 



1 Introduction 

It is well-known that a reliable quantitative prediction of the gluon fusion production cross 
section for Higgs bosons requires a next-to-next-to-leading order (NNLO) calculation (for 
a review on Higgs physics, see Refs. [H[2]). However, since it is a loop- induced process, its 
NNLO correction requires a three-loop calculation of a 2 — ► 1 process. Fortunately, it was 
found at next-to-leading order (NLO) [31 HJ [5] that the perturbative K-factor is very well 
reproduced in the so-called effective field theory (EFT) approach, where the gluon-Higgs 
coupling is taken into account by an effective Lagrangian 

£ c ff = —-^pxG^vG 1 ™ , (1) 

with H the Higgs field, G^ the gluonic field strength tensor, v = 246 GeV the vacuum 
expectation value of the Higgs field, and C\ a perturbatively evaluated Wilson coefficient 
(see, e.g., Ref. [U [7]). The NLO cross section in the EFT approach is then obtained by 
scaling the LO cross section (obtained in the full theory) with the effective NLO K-factor. 
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Although, to our knowledge, a quantitative understanding of the accuracy of this approach 
is still missing (in the sense that there is no error estimate) , the observed difference of less 
than 1% to the full NLO cross section (which is known in numerical form [8]) for Mjj < 2Mt 
was considered to be sufficiently convincing in order to trust the EFT approach also at 
NNLO. 

Apart from the inclusive NNLO calculation [9"1 110^ [TT]. the heavy-top limit has also been 
used for distributions, resummations, and even fully differential quantities at NNLO (for 
a review, see Ref. [T2] ). It is therefore of the utmost importance to justify the validity of 
the EFT approach. This has been first achieved at NNLO in Ref. [T31 Q31 Q21 Ql] by an 
expansion of the relevant Feynman diagrams in the limit My , s <C M t 2 , where y/J is the 
partonic center-of-mass energy. The apparent failure of this expansion for large s (which 
is only restricted by the hadronic center-of-mass energy squared, s) is only partly cured 
by the strong suppression of the parton luminosity. The prediction of the gg channel 
contribution, which accounts for more than 95% of the NLO hadronic cross section, was 
additionally treated by matching to the known large- s behaviour [T7]. The other channels, 
for which the large-s behaviour is not known, were treated by including only those terms 
in the 1/Mt expansion up to which the series was observed to converge. It was found that 
the resulting cross section agrees with the EFT result to better that 1% over the relevant 
mass range between 100 and 300 GeV. 

In this paper, we extend the analysis of Refs. [13 |, I15 |, ITT] by deriving the high-energy limits 
of the other channels as well. This leads to a significant stabilization of the qg channel 
which contributes about 2-5% to the total NLO cross section. The peculiar threshold 
behaviour of the quark-anti-quark channel prohibits a reasonable approximation from the 
high- and the low-energy information alone, but its contribution is in any case only at the 
per-mille level. We can therefore safely claim to present a stable prediction for the total 
cross section including top quark mass effects for Higgs masses between 100 and 300 GeV. 

2 Preliminaries 
2.1 Notation 

For the convenience of the reader, let us outline the notation at the very beginning. The 
Higgs mass is denoted by Mh, the on-shell top quark mass by Mt, and the hadronic and 
the partonic center of mass energies are s and s, respectively. Unless indicated otherwise, 
a s = (/^r) denotes the strong coupling in the MS scheme for five active flavours at the 
renormalization scale /xr. The following variables will turn out to be useful throughout 
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the text: 



Ml Ml 4M? s 

— • x = — ■ t= mK- " = 1 
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F = HI TT tR = ill TT = HI 7T , 

Ml H Ml Ml 
with the factorization scale //p. 

The inclusive hadronic cross section £jp P ' for Standard Modell Higgs production in proton- 
(anti-)proton collisions is obtained by convoluting the partonic cross section a a p for the 
scattering of parton a with parton (3 by the corresponding parton density functions (j> a / p (x) 
(PDFs): 

a a p(z,T,l F ) = / du}S a p(u,ti F )<T a p(z/u},T, l F ) , 



(3) 

c*,l3£{q,q,g} 



'(Z,T) = ^ <7al3(z,T,l F ) , P^{P,P}, 
c*,l3£{q,q,g} 

£ a p{u,m)= / [<t>a/p(y,HF)(f>f3/p'{u/y,flF)] ■ 

Juj y 



Note that the a a p depend on the factorization scheme (we use MS throughout this paper) ; 
only their sum a pp / is physical and thus formally independent of the factorization scale. 
Nevertheless, it will be useful to study the individual contributions to the total cross section 
separately as they have very different characteristics. The £ a g are parton luminosities and 
will be discussed in more detail in Section \'3.1[ 

We write the top quark induced partonic cross section as 

a a p(x, t, l F ) = cr (r) A a p(x, r, l F ) , (4) 

with 

(5) 



T^V^Gf (U s \ 2 2 , . 2 1 

1 + (1 — rj arcsin —= 



a (r)= " v r f ( — ) r 2 
256 V vr ' 



where Gf ~ 1.16637- 10 5 GeV 2 is Fermi's constant. The kinetic terms assume the form 



A a p(x,T,l F ) = S ag Sp g S(l - x) + ^2 ( — ) A^(x,T,ZF,fe) ■ 



TT 

n>l 



(6) 



At NLO, the full Mt dependence is known in numerical form [3] (the virtual terms are 
known analytically [TBI [T9l [20] ) . 
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Figure 1: Parton luminosities £(cj = z/x) at the LHC for Mt = 170.9 GeV at 
(a) M H = 130 GeV and (b) M H = 280 GeV, plotted as functions of x = M&/8. 
The vertical line denotes the threshold s = 4M t 2 . 



A fully general result for the partonic cross section at NNLO is as of yet unknown. In 
Refs. [131 f!4"l [T5l 116]. it was evaluated in terms of an expansion of the form 



/ M 2 \ 1 

A a p(x,T,fJ, F ) = ^2 [Y'l I) tt al3j i(x,lt,fJ,F) 



(7) 



with the analogous perturbative expansion as in Eq. ([6]). At NNLO, the first four terms 
(i < 3) have been evaluated \13\ I14j. The so-called EFT approach which has been used in 
all higher order analyses up to now, can be derived from the leading term of this expansion: 



(2, it) 



a,/3e{g,q,g}'' z 



(8) 



where an is given in Eq. ([5]). 
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Figure 2: Contributions of the partonic to the hadronic cross section from below 
(s < 4:M^; dashed) and above (s > 4M t 2 ; dotted) threshold, for the gg, qg, and 
the qq channel at NLO (gg includes the LO contribution). Note that qg uses a 
linear scale, while for the gg and the qq it is logarithmic. 
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3 Large- s limit 



3.1 Outline of the problem 

The expansion of Eq. ([7J) is expected to converge within s,M^ < 4M t 2 . While the Higgs 
mass range implied by electro-weak precision measurements lies comfortably in this range, 
the partonic center-of-mass energy y/J reaches values far beyond it, both at the LHC and 
the Tevatron. The corresponding breakdown of convergence manifests itself in inverse 
powers of x = M^/s, arising from 



M ? M t 



(9) 



Thus, in general, 



i ~ — as x — > . (10) 

x l 

Note, however, that at small x = M^/s there is a strong suppression by the parton 
luminosity S a g of Eq. ([3]) which we display for the various sub-channels in Fig.[TJ This, 
together with the fact that ^ a /3,o has no power singularities as x — > 0, are the main reasons 
that the heavy-top limit defined in Eq. ([8]) works so well. 

A further illustration of this observation is shown in Fig. [2] which compares the contribu- 
tions to the hadronic cross section arising from below (s < 4Af 2 ) and above threshold for 
the various subchannels at NLO. For the dominant gg channel, the region above threshold 
contributes only of the order of 2%. 

However, the spurious l/x singularities described before imply that in order to improve 
on the heavy-top limit by including higher terms in 1/M f , one needs to incorporate infor- 
mation on the large-s region. Fortunately, the leading terms can be obtained from general 
considerations. In the case of the dominant ^-channel, this was done in Ref. [T7j. This 
result was then combined with the 1/Mt expansion in Refs. |13^ I15j. 

Considering Fig.[H it appears that the center of the qg luminosity is at significantly lower 
values of x = M^/s than for gg. Correspondingly, the influence of the region above 
threshold is larger, as can also be seen in Fig.[2j The proper treatment of this region is 
thus much more relevant in the qg case. In addition, it is clear a priori that the EFT 
approach, which assumes that the top mass dependence at higher orders is determined by 
the LO one, cannot work as well in the qg channel which occurs only at NLO. In fact, the 
contribution of the qg channel to the total cross section in the EFT differs from the exact 
result by roughly a factor of two in the mass range between Mh = 100 and 300 GeV [8]. 
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In the next section, we extend the analysis of Ref. [T7J to the qg and pure quark channels 
(qq, qq, qq') at NLO and NNLO. The combination with the results of the 1/Mj expansion 
is done in Section [331 



3.2 Derivation of the leading high-energy behaviour 

The procedure to compute the leading logarithmic behaviour (LLx) of the partonic coef- 
ficient function to all orders in the strong coupling a s is based on /c^-factorization |21j . 
This technique has been used to resum coefficient functions for a few processes, e.g. heavy 
quark production [22], [23] , deep inelastic scattering [23] , Drell-Yan processes [25] and di- 
rect photon production [26J. The small x behaviour of Higgs production in gluon fusion 
was first computed in Ref. [27], in the heavy top approximation. The case of finite top 
mass was considered in Ref. |17j . where it was shown that, as expected, the coefficient 
function has only single high energy logarithms, while double logarithms appear in the 
effective theory. In Ref. j!7j , and in the phenomenological analysis of Ref. [28], only the 
gluon-gluon channel was considered. In the following the small x behaviour of all the other 
channels is computed, using high energy colour charge relations. For the sake of clarity, 
we set hf = //r throughout this derivation. 

The partonic cross section which enters the /c^-factorization formula is the leading order 
cross section for the process gg — > H, computed with two incoming off-shell gluons of 
momenta ki^, with kf 2 = ~ Iki^l 2 , contracted with eikonal polarizations. The impact 
factor is defined as the triple Mellin transform of the off-shell cross section 

h(N,r,M 1 ,M 2 ) = MM J^dCC"- 1 J" Wf 1_1 j™<%2^- X £ (11) 
MV ff (C,T,£i, 

where 



« - Tf2~' C - 7T77 Z TT T7T \ 12 ) 



Ml 

M' 2 H ' ' 2(h ■ k 2 - ki • k 2 ) 

and (p is the angle between the transverse polarization vectors ki and k 2 . In Mellin space 
the high energy limit corresponds to N — > 0; moreover it is easy to see that Mj — > is 
the collinear limit. The leading high energy behaviour of the coefficient function in the 
gluon-gluon channel is then found by identifying 

Mi = M 2 = ls , (13) 



N 

where j s is the anomalous dimension which is dual to the LO BFKL kernel xo> i- e - 

Xo (7s(a s /N)) = — , (14) 
a* 
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fc=i v 7 



(15) 



To all orders in perturbation theory, the leading logarithmic contribution to the MS coef- 
ficient function is 



Agg{N, T, /i F ) = h(0, T, 7s , 1S )R 2 {1S) 



M 



2 \ 2 7s 
H 



l4 



(16) 



The factor R is a scheme dependent function, first computed for MS in [23] ■ A recent 
calculation [29] has questioned that result. Although this issue must be solved for the 
resummation of the small x logarithms, it is not relevant for our present discussion. Our 
target is to compute the LLx behaviour of the coefficient function through NNLO, but the 
scheme dependence starts only one order higher: 



R = 1 + 



(17) 



The high energy behaviour of the other partonic channels can be derived from the gluon- 
gluon one by noticing that at LLx we have 



Cf 
C a 



Is , Jqq ~ Iqg ~ . 



(18) 



This means that, at LLx, a quark may turn into a gluon, but, because jgg is next-to-LLx, a 
gluon cannot turn into a quark. This leads to the following relations between the partonic 
coefficient functions and the gluonic impact factor [24] : 



A qg (N,r,fi F ) = -f / l (0,r, 7s ,7 s )i? 2 ( 7s ) 



M 



2 \ 27s 
H 



(19) 



C F 



f4 



A qq (N,T,fi F ) = -f h(0,T, ls , ls )R 2 ( T 



Ca 

2/i(0,r, 7s ,0)i?( 7s ) 



M 



2 \ 27s 
H 



i4 



(20) 



M 2 \ 7s 

-f +h(0,T,0,0) 



Notice that in the high energy limit A qq = A qq i 
identical and qq' to the distinct flavour case. 



*qq 



N- 



, where qq refers to the 
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The impact factor can be expanded in powers of Mj, which corresponds to an expansion 
in powers of a s : 



h(0,T,M u M 2 ) = M 2 H a Q {T)\l + h (1 \T){M l + M 2 ) 

+ /i (2) (r)(M 1 2 + Ml) + h^ 1) {r)M 1 M 2 + 



(21) 



The coefficients hP~\ and /i^ 1,1 ) have been evaluated numerically in [17^1 . The only 
difference here is that, in order to compute the LLx behaviour of all partonic subprocesses, 
we must keep the contributions M 2 ) and h^ 1 ' 1 ^ separated. 

It is then easy to substitute Eq. (j2"Tj) into Eqs. (fT9|) . (j20|) and invert the N Mellin transform 
to obtain the result in x space. Through NNLO, the small x limit of the partonic coefficient 
functions can be written as follows: 

A gg (x,T) = 5(l-x) + ^ \b$(t) - 2C A l F + O(x) 

+ (f) 2 [K? W - MaBV{t)I f + 2C\l F ) In \ + Bg)(r) + <9( a 



A M (x,r) = ^ fflW(T) - + 0(x) 

7T L ,3 



+ 



a. 



7T 



4? W - |cfp.B« (r)Z F + \c A C F l%\ In 1 + flg>(r) + 



(22) 



A gg( x ' r ) = A gg(^ r ) = A 99 '(a;,r) = 

2 r ^ )(r) " IH )(r)/F + ln l + ^ )(r) + ° {x) 



{Os_ 
V 7T 



Recall that we set fip = /ir in this section. The full fip, /UR-dependence is obtained by 
replacing 



«s(/4) 



7T 



— A)1fr + 



7T 



((A)^fr) 2 — /3Jfr) 



(23) 



in a A a/3 , where Z F R = ln(/4//4.) and /3 = 23/12, ft = 29/12. 

The coefficients A^l and are provided in the form of numerical tables in Table [TJ 
For all coefficients, the dependence on r is very smooth and can safely be interpolated 

(2) 

by straight lines, for example. The NNLO constants Bi are currently unknown; their 
influence on the final result will be studied at the end of Section HI 



1 See Eq. (36) and Eq. (38) of that paper, but notice the differences in the notation, e.g. the definition 



of r. 
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r 


n aa 


r>(l) 

D aa 




r 


A aa 


A ia 


Ai 2 ) 


1.0 


-0.8821 


-0.1960 




1.0 


33.0465 


8.7703 


1.2681 


1.5 


2.9212 


0.6492 




1.5 


35.9907 


9.5484 


1.3782 


2.0 


5.0234 


1.1163 




2.0 


44.2884 


12.2677 


2.1563 


2.5 


6.5538 


1.4564 




2.5 


53.1336 


15.1924 


3.0088 


3.0 


7.7650 


1.7255 




3.0 


61.8029 


18.0679 


3.8524 


3.5 


8.7693 


1.9487 




3.5 


70.1088 


20.8272 


4.6644 


4.0 


9.6279 


2.1395 




4.0 


78.0127 


23.4553 


5.4393 


4.5 


10.3781 


2.3062 




4.5 


85.5245 


25.9547 


6.1771 


5.0 


11.0444 


2.4543 




5.0 


92.6698 


28.3331 


6.8798 


5.5 


11.6437 


2.5875 




5.5 


99.4782 


30.6002 


7.5501 


6.0 


12.1883 


2.7085 




6.0 


105.9788 


32.7654 


8.1907 


6.5 


12.6875 


2.8194 




6.5 


112.1985 


34.8374 


8.8039 


7.0 


13.1482 


2.9218 




7.0 


118.1616 


36.8244 


9.3922 


7.5 


13.5760 


3.0169 




7.5 


123.8897 


38.7333 


9.9576 


8.0 


13.9752 


3.1056 




8.0 


129.4021 


40.5706 


10.5019 


8.5 


14.3495 


3.1888 




8.5 


134.7158 


42.3420 


11.0268 


9.0 


14.7018 


3.2671 




9.0 


139.8461 


44.0524 


11.5337 


9.5 


15.0345 


3.3410 




9.5 


144.8064 


45.7063 


12.0240 


10.0 


15.3497 


3.4110 




10.0 


149.6090 


47.3077 


12.4989 


10.5 


15.6491 


3.4776 




10.5 


154.2646 


48.8603 


12.9594 


11.0 


15.9343 


3.5410 




11.0 


158.7829 


50.3672 


13.4064 


11.5 


16.2065 


3.6015 




11.5 


163.1728 


51.8315 


13.8408 


12.0 


16.4670 


3.6593 




12.0 


167.4422 


53.2556 


14.2633 



Table 1: Coefficients for the large- s behaviour at NLO (left table) and NNLO (right). 
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3.3 Merging and partonic results 



Let us recall the knowledge of the partonic cross section at NNLO. Below threshold (s < 
4M t 2 ), the result is known in terms of an expansion in 1/M t and (1-x) [H Q22 [H@. Both 
expansions are expected to converge very well as long as Mh < 2Mt ~ 340 GeV. This is 
indeed observed for the gg and the qg channels at NLO in Fig. [3] for Mh = 130 GeV and in 
Fig. [H for Mh = 280 GeV. They show the partonic cross sections below threshold, keeping 
terms of order (1 — x) a (l/M^) b . In the left columns, a = 0, ...,8 and 6 = 5 (long to 
short dashes), while in the right columns, a = 8 and b = 0, . . . , 5. These figures compare 
the expansions to the exact result which we derived using standard techniques (see, e.g., 
Ref. [30J). In fact, the behaviour of the expansions suggests that, below threshold, the 
final result for the gg and the qg channels is numerically almost equivalent to the full Mt 
and x dependence. 

The NLO qq channel, on the other hand, has a very peculiar structure at threshold. At 
this order only one diagram with qq annihilating into an s-channel gluon contributes. 
Such a diagram is not enhanced in either the large- or small- x region, leaving room for 
a relatively pronounced structure at the threshold which cannot be described properly in 
our approach. However, the contribution of the qq channel to the hadronic cross section 
is down by almost three orders of magnitude relative to the gg channel, and still a factor 
of ten relative to the qg channel. We will nevertheless investigate its influence on the final 
prediction in more detail below. At higher orders we expect this effect to be reduced, 
because other diagrams with non-trivial high- or low-x limits will contribute. 

The corresponding curves at NNLO are shown in Figs.lSHHJ There is no exact result that 
one could compare to, but the quality of the convergence both of the 1/Mt and the (1 — x) 
expansions below threshold convincingly shows that they approximate the exact result to 
a very high degree in this region. 

From Section [3. 21 we know the leading high energy behaviour for general values of Mt and 
Mh- There are many ways then to merge the available information into a smooth function 
with the correct high- and low-energy behaviour (see, e.g., Refs. [T71 [281 02] ) • We decide 
to use [13] 



N 

k 



^fc- X 

k=l 



(24) 



+ (l-x)^ [aoBg-*W,( ) 

where cr^/p N^ x ) denotes the soft expansion of the partonic cross section through order 

2 Recently, the full x dependence was derived [16| . However, as argued before, the x dependence of the 
1/M t expansion does not hold for x < Mn/(4M t 2 ). 



11 



V°o (NLO) 




V°o ( NL0 > 



M H =130GeV 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

x 




M H =130GeV 



<W°o (NLO) 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

x 



Oqqb /0 O (NLO) 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

x 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

x 



<W°o ( NL0 > 



M H =130GeV 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

x 

V/°o (NLO) 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

x 



Figure 3: Partonic cross section at NLO for M t = 170.9 and Mu = 130 GeV 
at various orders in the expansion parameters (increasing order corresponds to 
decreasing dash size of the lines). Left column: C(l/M t 10 ) and 0((1 — x) n ), 
n = 0, .... 8. Right column: C((l - x) 8 ) and 0{1/M? n ), n = 0, . . . , 5. Solid: 
exact. The dashed vertical line indicate^the threshold. 
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Figure 4: Same as Fig.EH but for M H = 280 GeV. 
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Figure 5: Partonic cross section {gg and qg channel) at NNLO for Mt = 170.9 and 
Mh = 130 GeV for various orders in the expansion parameters (increasing order 
corresponds to decreasing dash size of the lines). Left column: C(l/M t 10 ) and 
0((l-x) n ), n = 0, . . . ,7. Right column: 0((l-x) 8 ) and C(l/M t 2n ), n = 0, . . . ,2. 
Solid: 0(1/Mf) and 0((1— x) 8 ). The dashed vertical line indicates the threshold. 
For the qq and the qq channel, see Fig.[6j 
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Figure 6: Same as Fig. 02 but for the qq and the qq channels (identical quark 
flavors). The figure for the qq' channel (different quark flavours) is not shown 
since it is almost indistinguishable from the one for qq. 
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Figure 7: Same as Fig. El but for M H = 280 GeV. 
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Figure 8: Same as Fig. El but for M H = 280 GeV. 
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A A 



<W°o (NLO) aja (NLO) 




Figure 9: NLO partonic cross sections as constructed from Eq. (|24h (with N = 8) 
by including successively higher orders in \/Mf . Dashed: 0(Mf n ), n = 0, . . . , 5. 
Solid: exact. Left/right column: M H = 130GeV/M H = 280 GeV. The dashed 
vertical line indicates the threshold. 
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a qg /a (NNLO) 




Figure 10: NNLO partonic cross sections (gg and qg channel) as constructed from 
Eq. (|24p (with N = 8) by including successively higher orders in \/Mf . Dashed: 
0(M t 2n ), n = 0,1,2. Solid: n = 3. Left/right column: M H = 130GeV/M H = 
280 GeV. The dashed vertical line indicates the threshold. For the qq and the qq 
channel, see Fig.QTJ 
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Figure 11: Same as Fig. llOl but for the qq and the qq channels (identical quark 
flavors). The figure for the qq' channel (different quark flavours) is not shown 
since it is almost indistinguishable from the one for qq. 
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(1— x) . For the unknown constants at NNLO we use the default values ctqB^ = a™n > ^(0), 
but we will study their influence on the NNLO hadronic results at the end of Section SJ 

At NLO, the resulting partonic cross sections are shown in Fig. O for Mh = 130 GeV and 
Mh = 280 GeV. The precision to which the gg and qg channels reproduce the exact result 
is quite impressive. As expected, the qq channel is approximated only very poorly though. 

The corresponding plots at NNLO are shown in Fig. [10] and \TT\ for the default values of 

(2) 

the high energy constant B^L Note that due to this undetermined constant, the curves 
do not all converge to the same point for x — > as is the case at NLO. Only the slope is 

(2) 

determined by the logarithmic coefficient AL. Nevertheless, convergence for the gg and 
the qg channel is very good, in particular at low Higgs masses. At larger Higgs masses, 
the x dependence becomes more and more unreliable as more terms in 1/Mj are included. 
This leads to the observed variations of the final result at Mh = 280 GeV. As will be 
shown at the end of Section 0] however, the variations affect the hadronic cross section by 
less than 1%. 

As expected, the qq channel does not seem to converge, but its contribution to the hadronic 
cross section is negligible as are those of the qq and the qq' channels (equal and differ- 
ent quark flavours, respectively). The observed convergence of the latter is much better 
though. In the following, we will include them in the total hadronic cross section, but we 
will only discuss the qq channel as representative of the pure quark channels (the one with 
the worst convergence behaviour). 



4 Hadronic Results 

In order to study the effect of the 1/Mt terms on the hadronic cross section, we define (see 
also Ref. [13]) 



a^°(Mn = a 6 ag 6p g 6(l - x) + (Mf ) 



a™ LO (Mf) = a S ag 6 Pg S(l - x) + + a™{M? 



(i) 



(25) 



where <t^(M™) is the N fc LO contribution to the partonic cross section evaluated as an 
expansion through (D(l/M™), and matched to the \ow-x limit as described in Section [3.21 
A^g ^ is the EFT result as defined in Eq. ([8]). Note that this differs from an extended EFT 

approach, where A^j would be expanded in terms of 1/Mt, while the full r dependence 
in <to(t) is kept. We will return to this latter approach at the end of this section. The 
corresponding hadronic quantities derived from Eq. (|25|) are denoted by <T^g°(M t n ) and 

a™ LO (M t «). 



21 



Fig.ll2lshows the relative gg, qg and qq contribution crjjg (M") to the total hadronic cross 
section. The dashed lines correspond to successively higher orders in 1/Mt, while the solid 
line shows the exact result. The curves are all normalized to the exact NLO cross section 
0.NLO p or gg and the qg channels, one observes excellent convergence towards the 
exact result (solid line) . The small deviations are reflections of the deviations between the 
solid and the dashed lines in Fig.[9](a)-(d). As pointed out above, the mass effects in the qg 
channel are quite large, ranging from roughly a factor of two to four in the relevant Higgs 
mass range. Of course, the overall size of the qg channel is below 5%. As expected, the 
picture in the qq channel is significantly worse. Only the order of magnitude is captured, 
but there is no sign of convergence towards the exact result whatsoever. Its contribution 
to the total cross section is only of order 10 -3 though and thus irrelevant. 

The corresponding plots at NNLO are shown in Fig.[T3J Since there is no exact result in 
this case, we normalize the curves to the full NNLO EFT result, cf. Eq. ([8]). Also, the solid 
lines always refer to the subchannels evaluated in the EFT approach. The observations are 
quite similar as at NLO: the difference between the EFT result and the 1/Mt expansion 
for the gg channel is about 1% which is of the order of the accuracy to which we expect 
the capture the mass effects. In the qg channel, the relative difference between the EFT 
result and the 1/Mt expansion is significantly larger (~ 20%), but the influence of this 
effect on the total cross section is again only of order 1% due to the strong suppression of 
the qg channel. The qq channel does not seem to converge very well, but is numerically 
negligible (the true mass effects are not expected to change this). 

The hadronic results for the Tevatron are shown in Fig. [TH at NLO and NNLO. The conclu- 
sions are very similar to those for the LHC, thus justifying the use of the EFT approximation 
for Higgs searches also in this case |31j . 

Overall, we conclude that the final result for the NNLO cross section including top mass 
effects is within 1% of the EFT result. 



(2) (2) 

Dependence on B a ^. As pointed out above, the constants B a ^ for the large- s be- 
haviour are currently unknown. From the curves in Fig. [10] and [11] our choice oqB^ = 
cr^g(O) seems to be reasonable, leading to rather smooth curves over the full x-range. 



Nevertheless, in order to estimate the uncertainty induced by this unknown constant, we 
set (JoBgg = t x <jgg (0) and find that the dependence of the hadronic cross section on t is 
very well described by a linear function: 

a NNLO „ (1 _ Q 01t) ^NNLO (2g) 
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Figure 12: (a)-(c) Sub-channel contributions to the hadronic cross section at NLO, 
normalized to the full result. Note that gg includes the exact LO contribution, 
cf. Eq. (|25[) . Dashed: including terms of order 1/Mt 2n in the numerator (n = 
0, . . . , 5 from long to short dashes). Solid: exact, (d) Sum over all sub-channels. 
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Figure 13: (a)-(c) Sub-channel contributions to the hadronic cross section 
at NNLO, normalized to the full NNLO EFT result (LHC conditions). Note 
that all channels include their lower order contributions in the EFT approach 
(cf. Eq. (|25p ). Dashed: including terms of order 1/Mj 2n in the numerator 
(re = 0,1,2,3 from long to short dashes). Solid: EFT result, (d) Sum over 
all sub-channels. 
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Figure 14: Hadronic cross section at (a) NLO and (b) NNLO for the Tevatron, 
normalized to the full NNLO EFT result. Note that the lower order contributions 
are included in the EFT approach (cf. Eq. (|25p ). Dashed: including terms of order 
1/Mt 2n in the numerator (n = 0, 1,2,3 from long to short dashes). Solid: EFT 
result. 



Again recalling the smoothness of the curves in Fig. [10] and [TT| we do not expect the 
parameter t to be significantly larger than one. The resulting uncertainty is therefore at 
most at the percent level and therefore much smaller than the scale uncertainty of the 
NNLO result. 



Is the heavy-top limit a coincidence? Let us conclude this section with a remark 
on the extended EFT approach as mentioned in the discussion after Eq. (|25p . It would 
be possible that the high quality of the EFT approach is a coincidence, in the sense that 
there is an accidental cancellation among the higher order terms in the 1/Mt expansion of 
the A a p. This would have a significant effect on the applicability of the EFT approach to 
other quantities, of course. 

However, we have checked that this is not the case. All the curves of the extended EFT 
approach lie within 1% of our final result. 
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5 Conclusions 



The hadronic Higgs production cross section due to gluon fusion was presented including 
effects from a finite top quark mass. We have extended previous analyses by deriving 
the high-energy limits of all partonic sub-channels and combining them with the known 
1/Mt expansions. Although the mass effects on the absolute size of the qg channel are 
large, they have no significant effect on the total hadronic cross section. Therefore, the 
main conclusions of previous analyses |14^ [16] remain valid, and the EFT approach is still 
justified. 
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